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Abstract 

We study a model in which a Hubbard Hamiltonian is 
coupled to the dispersive phonons in a classical non- 
linear lattice. Our calculations are restricted to the 
case where we have only two quasi-particles of oppo- 
site spins, and we investigate the dynamics when the 
second quasi-particle is added to a state correspond- 
ing to a minimal energy single quasi-particle state. 
Depending on the parameter values, we find a num- 
ber of interesting regimes. In many of these, discrete 
breathers (DBs) play a prominent role with a local- 
ized lattice mode coupled to the quasiparticles. Sim- 
ulations with a purely harmonic lattice show much 
weaker localization effects. Our results support the 
possibility that DBs are important in HTSC. 
PACS: 71.38.-k, 63.20.Pw and 63.20.Ry 



cumulating experimental evidence exists for electron- 
lattice effects in high temperature superconductors 
1101 ITT] , and theories based on electron-phonon in- 
teractions have also been proposed |21E1EJE]. Here 
we follow the idea that both electronic correlations 
and electron phonon interactions are important |12) 
and study a model in which a Hubbard Hamiltonian 
is coupled to dispersive phonons. Our main aim is 
to explore one extra ingredient, which has generally 
been ignored until now, the importance of the an- 
harmonic character of lattice vibrations. Whilst our 
ultimate aim is to understand HTSC, here we propose 
a specific mechanism for pair formation that involves 
the interaction of polarons through a nonlinear lat- 
tice mode, which will have applications in other areas. 
We study the stability of such a pair as a function of 
the electron-electron (or hole-hole) interaction. 



1 Introduction. 

In spite of the many studies Ql El El El E] made since 
it was first discovered high temperature super- 
conductivity (HTSC) remains a challenge. The na- 
ture of the carriers and the mechanism behind pair 
formation are still unclear. According to many re- 
searchers, HTSC can be explained by a purely elec- 
tronic model, such as that described by the t — J or 
the Hubbard Hamiltonians, for which charge and/or 
spin interactions are paramount. This view is essen- 
tially based on the absence of isotope effects seen in 
some experiments and the apparent ci-symmetry 
of the superconducting wavefunction. However, ac- 



2 The Hubbard-Davydov 
Hamiltonian. 

The Hamiltonian H we use has three parts: 

H = Hqp + ff qp _ pll + H ph ( 1 ) 

where -Hqp is the Hamiltonian for a quasiparticle 
with spin i, ^qp_p n describes the interaction of the 
quasiparticle with the lattice and is the lattice 
(phonon) Hamiltonian. 

The Hamiltonian for the quasiparticle is the ID 
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Hubbard Hamiltonian: 



The phonon Hamiltonian is as follows: 



n.cr n 



where the sums are over the sites n, going from 1 to 
AT, (N is the total number of lattice sites) and a refers 
to the spin and can be up or down. $ na is the cre- 
ation operator for a quasiparticle of spin a at site n. 
e is the self-energy of the quasiparticle, t the transfer 
term for the quasiparticle to move between neigh- 
bouring sites. We depart from the usual notation in 
that the on-site quasiparticle-quasiparticle coupling 
is here designated as 7 (and not U) to avoid confusion 
with the variables {u n } used for lattice displacement 
(see below). Both negative and positive values of 7 
will be considered, corresponding to the attractive 
and repulsive Hubbard models, respectively. 

As in the Davydov model for energy transfer in 
proteins ^qp-ph' ^ ne Hamiltonian for the inter- 
action of the quasiparticle with the lattice includes 
the coupling to acoustic (or Debye) phonons: 



H, 



qp-ph 



where x is a parameter which describes the strength 
of the quasiparticle-lattice interaction. Many previ- 
ous publications have included electron-phonon inter- 
actions in the framework of the model of Holstein ^1] , 
in which only short-range interactions are are consid- 
ered. As has been pointed out elsewhere ^5], when 
the electron screening is poor, such as in cuprates, 
electron-phonon interactions are long range, which 
can be described by acoustic phonons. 
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where u n is the displacement from equilibrium posi- 
tion of site n, p n is the momentum of site n, a is the 
equilibrium distance between sites, n is the elasticity 
of the nonlinear lattice and k' is a similar constant 
for the on-site potential. Here, the coupling interac- 
tions between sites are described by a Lennard- Jones 
potential, a potential commonly used to de- 

scribe interactions between atoms. In a high tem- 
perature cuprate, this potential describes the inter- 
actions of the copper and oxygen atoms in one Cu-0 
layer. The on-site potential is as used in many 

breather studies ^Hl- It can be considered to repre- 
sent the effect, in a mean field approach, of the rest of 
the crystal on the one dimensional chain whose states 
are studied explicitly. In a cuprate, this models the 
effect of the neighbouring layers on the Cu-0 layer. 

Our Hamiltonian (111 111 includes two sources of non- 
linear effects. The first comes from the intrinsic non- 
linearity of the Lennard- Jones potential, and the 

on-site potential, H®^. The second source of nonlin- 

earity is extrinsic and comes from the interaction of 
the quasiparticle with the lattice (cf. Eq. OJ). The 
former is the source of nonlincarity in the studies of 
discrete breathers JB] and the latter is the cause of 
localization in polaron theory. 

We adopt a mixed quantum-classical approach in 
which the lattice is treated classically, while the 
quasiparticle is treated quantum mechanically. Ac- 
cordingly, the displacements u n and momenta p n are 
real variables. The quasiparticle variables are opera- 
tors, a distinction which is marked by the hats above 
the operators. The importance of quantum effects 
of the lattice can be assessed by considering the full 
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quantum model at finite temperature, which has al- 
ready been done for the Davydov Hamiltonian. It 
was found that, at 0.7 K, the lattice displacement 
correlated with the position of the quantum particle 
in exact semiclassical Monte Carlo simulations dif- 
fered by 15 % from the corresponding variable in 
exact simulations in the fully quantum system. At 
11.2 K, the two approximations lead to virtually the 
same value J7|- We would like to emphasize that 
the approximation we consider here is not an adia- 
batic approximation. In an adiabatic approximation 
the kinetic energy of the phonons is neglected with 
respect to the kinetic energy of the quantum particle. 
We do not do that here, as our dynamical equations, 
eqs (6) and (7) below, include the time derivative 
of the momenta of the lattice sites. What we do is 
to consider that the dynamics of the lattice can be 
treated classically. Both the semiclassical (or quan- 
tum/classical, as we prefer to call it to differentiate 
from other use in the literature) approach we apply 
here and the adiabatic approximation lead to similar 
results when we consider the ground states of the sys- 
tem (because the corresponding solutions have zero 
kinetic energy), but they are different when we deal 
with dynamics, as we do in this work. 

Ultimately, the need for a full quantum treatment 
comes from comparison with experimental results. 
Isotopic effects can only be described in a fully quan- 
tum framework. Our main aim here is to explore the 
importance of anharmonicity in the lattice for the 
dynamics of paired states, something which is much 
more complicated to do within a fully quantum for- 
malism. Thus, as a first approximation, we restrict 
ourselves to the mixed quantum-classical regime and 
study the behaviour of a pair of quasiparticles, cou- 
pled to a nonlinear lattice. 

With these assumptions, the exact two quasiparti- 
cle wavefunction for the Hamiltonian (11IH is: 

\m)= E <WM,K},*)4 T 4iio> ( 5 ) 

n,m— 1,N 

where <f> nm is the probability amplitude for a quasi- 
particle with spin up to be at site n and a quasipar- 
ticle with spin down to be at site m. The proba- 
bility amplitude is dependent on the lattice displace- 
ments and momenta in a way that is not specified 



a priori and is determined by the equations of mo- 
tion. Similarly to other systems [T^], the equations 
of motion for probability amplitudes <j) nm are derived 
by inserting the wavefunction 10) in the Schrodinger 
equation for the Hamiltonian (121 II) . and the equations 
for the displacements and momenta are derived from 
the Hamilton equations for the classical functional 
£ 2 = {ip\H\il>}. They are: 




= -* + 4>3+U + <Pji-i + <frji+i) + 7<Pji$jl + 



+ x {uj+l - Uj-i + ui+i - ui-i) 4>ji (6) 
dt duj 

-x(\?U\ 2 -\v} +1 \ 2 + \vU 2 -\vU 2 ) 

where | ipj | 2 , the probability for the quasiparticle with 

spin up to be in site j and I'/'jl 2 , the probability for 
the quasiparticle with spin down to be in the same 
site. These are given by: 

N 

\<Pj I 2 = McjrCjiM = E l^'l 2 ' 

2 = 1 
N 

1=1 

3 Dynamical states. 

We consider the case in which the quasiparticle den- 
sity is low and the starting point is that of an iso- 
lated quasiparticle interacting with the lattice. We 
wish to find if the addition of a second quasiparticle 
with opposite spin to that state can lead to pairing 
of the two quasiparticles, and how the relative stabil- 
ity of the paired state depends on the quasiparticle- 
quasiparticle interaction 7. 

We start from the state of a single quasiparticle. 
The wavefunction is 

|#>=5>n fi n|J°> (8) 
n 

Minimum energy states for this one quasiparticle can 
be found by numerical minimization of the energy 



3 



functional E 1 — IV' 1 } with respect to the prob- 

ability amplitude for a single quasiparticle in site n, 
(j>\, and to the displacements u n [TJJ. Two kinds 
of minimum energy states are found. For sufficiently 
large quasiparticle-lattice interaction x, the quasipar- 
ticle states are localized and there is an associated 
lattice distortion. We call this the single particle po- 
laron, or simply polaron. Below a threshold value for 
X, the states are delocalized, as in the usual Bloch 
states, and the lattice is undistorted. We have con- 
sidered a value of % an d other parameters such that 
the initial one quasiparticle polaron state is neither 
too weak nor too stable when compared with delo- 
calized, Bloch states for the same values. While it is 
important to find the behaviour of the two quasipar- 
ticle states considered here for different values of the 
parameters, our choice ensures that the results here 
are not the consequence of extreme values. 

The dynamical states we study are the perturba- 
tions of the single polaron state, induced by the pres- 
ence of a second quasiparticle with opposite spin. Be- 
cause the number of variables (f> nm that characterize 
the wavefunction J5J increases with the square of the 
lattice size, in order to be able to integrate the equa- 
tions of motion for a sufficiently long time, the size of 
the lattice was kept relatively short, i. e. the number 
of sites is N = 20. The aim is to investigate the in- 
fluence of the strength and sign of the quasiparticle- 
quasiparticle interaction 7 on the dynamics of the 
paired quasiparticle states. 

The parameters of the simulations in the figures are 
the same, except for the quasiparticle-quasiparticle 
interaction 7. In Fig. ^ we se t 7 A = — 10 in an at- 
tractive Hubbard model. The addition of a second 
quasiparticle leads to a localized state for the pair, 
with a very slight peak oscillation, that is hardly vis- 
ible in the figure. (The probability for the second 
quasiparticle is the same as that shown and is not 
displayed). The lattice, however, sets into a breather- 
like oscillation |16| . i.e., a localized excitation with an 
internal oscillation. Indeed, at the site of the initial 
lattice distortion, oscillations are clearly visible in the 
lattice displacements and momenta. A striking ob- 
servation is that the amount of radiation generated 
is very small, and most of the energy of the lattice is 
associated with the breather. Fig. which displays 




Figure 1: Time dependence for (a) the probability 
for one quasiparticle to be in site n, (n = 1 • • ■ N, 
N — 20), (b) the lattice displacement and (c) the 
momentum of site n. Time is in picoseconds. The 
parameters are t — 10 x 10~ 22 J, \ — lOOpN, k = 
lN/m, k' = 2k, a = 4.5l and 7 = -100 x 10" 22 J. 
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another 6 ps period of the dynamics at a later time, 
demonstrates the stability of this solution. 

A Hubbard Hamiltonian with a much weaker at- 
traction, corresponding to a ratio of j/t = —0.5, is 
considered in Fig. |3 where the last 6 picoseconds of 
a 42 picosecond simulation are displayed. A mod- 
ulation of the peak of the probability distribution is 
now clearly seen, which has the same frequency as the 
main modulation of the lattice breather. The mod- 
ulation of the quasiparticle probability is associated 
with a periodic change of shape in which a lower peak 
with a slight tail appears. Even at this comparatively 
much weaker interaction, the amount of radiation is 
very small and most of the lattice energy is in the 
breather. The frequency of the main modulation of 
the breather is as for 7/t = —10. 

In Fig. 0| a repulsive Hubbard Hamiltonian is con- 
sidered, with 7/t = +0.5. The modulations and the 
associated tails of the probability distribution for the 
quasiparticle are now more pronounced, but their 
main frequency is unchanged. Although there is a 
slight increase in the radiation in the lattice, the sta- 
bility of the breather and of the quasiparticle solution 
is apparent. 

In Fig. the repulsive interaction is increased to 
7/t = +1. The modulations in the probability distri- 
bution for the quasiparticles lead to greater periodic 
changes of shape, still with the same frequency as for 
the other values of 7. The radiation in the lattice is 
now more visible, but the breather remains stable. 

In Figures IH1 and [7| a large repulsive value, corre- 
sponding to 7/t = 5 is taken. This leads to a change 
in the probability distribution for the quasiparticles, 
from a single site peak into a two site peak, with 
periodic oscillations which make one probability at 
one site larger than the other. The lattice variables 
show that, concurrently with the appearance of the 
breather, a considerable amount of radiation is gen- 
erated. Also noticeable is the fact that the frequency 
of the modulations has changed. Fig. [7| shows that 
the new quasiparticle probability distribution is sta- 
ble, as well as the lattice breather, even if the noise 
which results from successive passes of the radiation 
through the periodic boundaries, constitutes a signif- 
icant part of the lattice energy. 

In Figures |H1 and 03 a repulsive interaction corre- 
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Figure 7: Same as Fig. but at a later time. 



sponding to j/t — 10 is used. Fig. |S| shows that a 
drastic transformation takes place in which the ini- 
tial distribution changes into a two peak distribu- 
tion. One of the peaks is located where the initial 
lattice distortion was and the second peak is as far 
away from it as it can be in this lattice. Also, while 
the peak that is located at the original lattice dis- 
tortion site remains unmodulated in time, as well 
as its associated lattice distortion, the second peak 
oscillates with approximately the same frequency as 
that in Figs. Eland [3 The momenta in Fig. [S] show 
3fjearly that the second peak has an associated lattice 
breather, while the first peak is associated with a dis- 
tortion that is essentially static. After some time, be- 
cause of the repeated reflection of the radiation from 
the boundaries, this picture is not so clear. Both 
peaks show oscillations in the displacements and the 
momenta of the lattice are rather noisy. However, 
Fig. [5| does illustrate the stability of the two peak 
solution, even in the presence of such relatively large 
amount of noise. 

^% Dynamical states in the fully 
harmonic approximation 

The early theory of pair formation via interaction 
with phonons assumed that the lattice motion was 
harmonic. It is interesting to see how the dynamics 
of the two electron states would be in this case, and 
this section is devoted to that question. The first two 
terms (JSJ, © in the Hamiltonian we consider in this 
section are the same as before, but now the phonon 
3Hamiltonian is given by: 

^harm _ ^co-harm _|_ ^os-harm ^ 



ph JJ ph 
i N 

#co-harm = l K £ K _ Un _ i)S 

n=l 
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The phonon Hamiltonian @ can be obtained from 
(QJ by considering the limit of small displacements, 
in which only the linear terms remain. In this case, 
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the only nonlinear term left for the total Hamiltonian 
is that which describes the quasiparticle-lattice inter- 
action. It should be pointed out that, if we disregard 
the correlation term in the equations of motion 
for this system are those studied by a number of au- 
thors ^1 for a single single polaron, and for any 
number of polarons by Alexandrov |21j . 

Fig. 1 1UI shows that when the effective interaction is 
such that j/t — —10, the addition of an extra elec- 
tron to the minimum energy single polaron leads to a 
state in which both electrons are in the same site with 
a strong lattice deformation of breather type associ- 
ated with their presence. The time evolution of the 
momenta, however, shows that there is no breather 
formation, only phonons which travel along the lat- 
tice. Because of the periodic boundary conditions, 
these phonons eventually come back and after they 
have crossed each other many times the lattice be- 
comes very noisy. The lattice deformation associated 
with the two electrons oscillates periodically because 
of the interference of these phonons, but does not 
move. Also, the state of the two electrons remains 
localized on one site all the time. 

Similar dynamics takes place for 7/f = —0.5, ex- 
cept that very slight oscillations in the probability 
distribution for the electron states also takes place 
(not shown). 

When the electron-electron interaction is repulsive 
and such that 7/i = +5, the phonon emission leads 
to fluctuations in the electron probability distribu- 
tion that are clearly visible in Fig. ^2 The dynamics 
is similar to that of fig. 1101 with phonons propagating 
along the lattice and causing oscillations in the other- 
wise constant distortion induced by the two electrons. 
Again, the momenta show that there is no breather 
formation and all the dynamics of the lattice is due 
to the phonon propagation and interference. 

For a repulsive interaction for which j/t = +10, 
the two electrons split up and the probability dis- 
tribution shows two peaks, both of which have an 
associated lattice deformation with the breather pro- 
file (see Fig. I12t . Phonons are generated from each of 
these locations and their interference eventually leads 
to a noisy lattice. The two peaks in the probability 
distribution for the electrons oscillate in a less regular 
fashion than in the anharmonic lattice, but remain 




Figure 10: Same as Fig. ^ but with 7 = —100 x 
10~ 22 J and for the harmonic lattice 
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stable throughout the simulation. It should be no- 
ticed that for this harmonic approximation also, the 
lattice displacements induced by the electrons/holes 
are not small. Hence, an accurate representation of 
the dynamics should include the nonlinear terms in 
the lattice Hamiltonian, as was done in the previous 
section. 

5 Discussion 

Our aim was to investigate the relative stability of 
a correlated pair of quantum quasiparticles with op- 
posite spins with respect to their uncorrelated states. 
The starting point was thus the state of a single quasi- 
particle polaron and we studied the dynamic states 
which arise when a second quasiparticle is added to 
the first state. The Hamiltonian used includes several 
physical ingredients. On the one hand, it contains 
two sources of nonlincarity, one intrinsic to the lat- 
tice and another which arises from the quasiparticle 
lattice interaction. Such nonlinear lattices have been 
shown to possess generic solutions known as discrete 
breathers (DBs) ^5] ■ The study of systems in which 
nonlinear lattices are coupled to one quantum quasi- 
particle, on the other hand, is just beginning |22M23| . 

To our knowledge, this is the first time that the 
coupling of two quantum quasiparticles to a nonlin- 
ear lattice has been considered. Indeed, a second in- 
gredient is the inclusion of quasiparticle-quasiparticle 
interactions, in addition to the quasiparticle-lattice 
interactions found in the polaron model. The 
quasiparticle-quasiparticle interactions can represent 
Coulomb interactions, and/or spin-spin interactions, 
and be either attractive or repulsive. The dynamical 
simulations indicate that for these extended systems, 
DBs are generic solutions also and can be generated 
by the presence of a second quasiparticle. These lat- 
tice breathers can in turn stabilise localized, paired, 
quasiparticle states, for a large range of 7 values. 
Windows of 7 were found for which similar solutions 
are obtained. Thus, for a ratio of j/t between —10 
and +1 (Figs. IT15|l . DBs are found in the lattice and 
in the quasiparticle, with the same main modulation 
frequencies. For larger values of j/t, two different so- 
lutions were found (see Figs. l6l9|l . In one solution the 



quasiparticles distribution is split into equal values in 
two neighbouring sites and in the second a two peak 
distribution, with the peaks as far apart as possible 
in the lattice used, is observed. 

This Hamiltonian includes the two main physical 
causes for quasiparticle pairing that have been con- 
sidered in HTSC and allows for interpolation between 
them, by varying the strength of the relevant param- 
eters. According to our results, a greater importance 
of quasiparticle-lattice interactions in pair formation 
should arise in systems for which the dynamics of 
the lattice is fast enough compared to the quasipar- 
ticle dynamics, so that the lattice relaxes when the 
two quasiparticles meet. Conversely, a corresponding 
greater importance of quasiparticle-quasiparticle in- 
teractions should be associated with systems in which 
the lattice dynamics is much slower than the quasi- 
particle dynamics. 

An implicit assumption in this study is that the 
nonlinear character of the lattice plays an important 
role in HTSC. Although the lattice distortions are 
weak in conventional superconductors, and thus the 
lattice dynamics can be approximately described by 
a linear system, we argue that in HTSC these dis- 
tortions are such that the lattice enters a nonlinear 
regime. This may be why the sound velocity de- 
creases by a few parts per million in conventional 
superconductors, whereas in a high Tc material there 
is an increase which is two or three orders of magni- 
tude larger than in the former case. Our simulations 
with the harmonic lattice show that the percentage 
of energy transferred to travelling phonons is much 
larger than for the anharmonic lattice. 

The breather-like solutions found in the dynamical 
simulations are a signature of the nonlinear dynam- 
ics of the lattice. The possibility that breathers are 
associated with HTSC has been suggested elsewhere 
PIE]. Our study indicates that DBs are generic 
excitations in systems governed by the Hamiltonian 
used here. Moreover, within a certain range of the pa- 
rameters, the states in which two quasiparticles are 
paired and coupled to a DB are energetically more 
favourable than those of uncorrelated quasiparticles. 
Hence, this study gives weight to the possibility that 
DBs are important in HTSC. 
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